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(57) Abstract 

A system for locating objects in space, such as medical instruments (16) within the body of a patient (P), based upon transmission of 
magnetic fields from coils (10) in a fixed frame of reference to sensors (20, 22, 24) on the object or vice-versa. The current supplied to the 
coils is adjusted to assure that the sensors receive fields within a preselected range of magnitudes regardless of the location of the object 
in space. This assures that the sensor operates within its optimum range, and permits use of compact transmitters and sensors. 
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DESCRIPTION 

MAGNETIC LOCATION SYSTEM 
WITH ADAPTIVE FEEDBACK CONTROL 

TECHNICAL FIELD 

The present invention relates to systems for 
determining the location and/ or orientation of objects 
in space by detecting magnetic fields. 
BACKGROUND ART 

Various systems have been proposed for 
detecting the position and/ or orientation of an object 
using magnetic or electromagnetic fields. These systems 
15 typically employ field transmitters, such as 
electromagnet coils, disposed at known locations in a 
fixed reference frame and a sensor, such as a coil or 
other transducer mounted to the object to be located. 
Each transmitter projects a field varying in space in a 
20 fixed frame of reference. The pattern of variation in 
space for each transmitter is different than the pattern 
for each other transmitter. For example the 

transmitters may be identical to one another but 
disposed at different locations or in different 
25 orientations. The field patterns of the transmitters 
are thus displaced or rotated relative to one another 
and relative to the fixed frame of reference. The sensor 
on the object detects the parameters of the field 
prevailing at the location of the object as, for 
example, the magnitude and/ or direction of the field at 
the object or the magnitudes of individual components of 
the field at the object in one or more preselected 
directions. The transmitters may be actuated in a 
predetermined sequence so that at any time only one 
transmitter is active and therefore the field prevailing 
at the object is only the field contributed by one 
transmitter, plus a background field due to the Earth's 
magnetic field and other environmental sources. 
Alternatively, the transmitters can be driven at 
different frequencies so that components of the signal 
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from the sensor varying at different frequencies 
represent contributions to the field at the object fro- 
different transmitters. Based upon the detected 
parameters of the fields from ^ indiviaual 
transmitters, and the known pattern of variation of the 
field from each transmitter, a computer system 
calculates the position and orientation of the sensor, 
and hence the position of the object bearing the sensor, 
in the fixed frame of reference of the transmitters. In 

carrier! 0 ' " hlS *" *" 

blurV. , « transmitters, whereas a 

plurality of sensors are disposed at various locations 
or orientations in the fixed frame of reference. The 
location and/or orientation of the object is deduced 
from signals representing the parameters of the field 
prevailing at the various sensors. 

m n •* T tBmS ° f thlS general nature *** disclosed 
in united states Patents 4,849,692; 4,642,786- 
4,710 708; 4,613,866 and 4,945,305. Systems according 
to this general design can be used to provide a three! 
dimensional spatial input capability for a computer. 
Another system of this nature is disclosed in 
international Patent publication WO94/04938. In the 
-938 publication, the object to be located may be a 
medical endoscope. Such a system may include a sensor 
mounted on the tip of an endoscope, so that the location 
and/or orientation of the endoscope tip can be 
determined while the sensor is disposed inside the body 

30 mLt "T^ PatiSnt - aU ° WS *»» Physician to 

monitor the endoscopic procedure without resorting to a 
fluoroscopy or other techniques using ionizing radiation 

26 2 ,T , the inStrU * ent - In «• Pediment (pages 
26-27, the .938 publication contemplates an arrangement 

sL,? f transmltti "9 coils are actuated 

simultaneously to "steer., the direction o, the resulting 
field and thereby align the field with the sensor. The 
938 publication characterizes this arrangement as 
undesirable. other systems for locating medical 
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instruments such as endoscopes and catheters based upon 
transmitted fields are disclosed in United States Patent 
No. 5, 042,486; 5,099,845; 5,211,165; 5,251,635; 
5,253,647; 5,255,680; 5,265,610 and 5,391,199. 

Systems of this nature typically have used 
fields having strength which varies as the third power, 
or a higher power, of distance from the transmitter. 
There is, accordingly, a very large variation in the 
strength of each field from location to location. In a 
medical device locating system, adapted to detect the 
location of the sensor anywhere within a sensing volume 
having dimensions on the order of 0.5 to 1.0 m, the 
field strength may vary through many orders of magnitude 
from one end of the sensing volume to another. When the 
sensor and the object to be located happen to be near a 
particular transmitter, they will receive an extremely 
strong field, whereas when the sensor and object are at 
the end of the sensing volume remote from the particular 
transmitter, they will receive only a very weak field 
from that transmitter. This requires a sensor with a 
very substantial dynamic range, capable of accurately 
monitoring very strong and very weak fields. That, in 
turn, poses stringent demands on the design of the 
sensor and may require a larger sensor. Moreover, the 
signal to noise ratio is poor when monitoring a very 
weak field and the accuracy of the system is accordingly 
less than optimum in these conditions. Both of these 
considerations are particularly important in systems for 
medical applications, where the sensor must be small to 
fit inside a medical device such as a catheter or 
endoscope. Although such systems can provide useful 
results, further improvement would be desirable. The 
same problems arise in systems using a coil or other 
transmitter on the object to be located and plural 
sensors in the fixed frame of reference. Here again, a 
particular sensor will be exposed either to a very 
strong field or a very weak field depending on whether 
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Commonly assigned, copending United states 
Patent 0 S/132 ,47 9 discloses , sol ution to 

l: c r rtain prefer ~ d — ac^r 

operate J " lnClUae ^"^s *** — 

"eTds Bacrr?; 3 P1Urality ° £ 

J - Leias « Each field has a+- 

nonzero magnitude that Ts 8 

proves less vanatxon ln tl . w magnltude ^ 

w L ;:iir7 than a 

tho 016 <J»asl-linear fields, 

the difference between minimum a nd naxlmm ; 

• given sensing volume is substantially seller TWs 
substantially alleviates the problem discussed' ablve 
However, it is not always convenient to use The 

iS^ 1 " provide 
linearly varying fields. 4 

stoimmv n T tub Tflv E n .T. T 

one aspect of the present invention provides 
apparatus for determining position including rield 

":r ing r ans . for producing * pi< ™ « u:: : 

fields. The Held generating means is arranged to 
prov.de such fields having parameters, such as Ltd 
strength and/or field direction, varying with location 

variat.cn. The pattern of variation for each field is 

. "des r P i™ e ; :rr the £isia ~ tin * 

are disposed" a t ^ZZ^Z^ «" ^ 

« • ^ . airrerent locations and/or different 

35 -stations relative to the sensing voiume. Apparatus 
a coring to this aspect of the invention also Llu^ 
at least one sensor adapted to detect one or more 
parameters of the field prevail™ ... .... 

prevailing at the sensor when 
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the sensor is at an unknown location within the sensing 
volume, and to provide one or more sensor signals 
representing such detected parameters. For example, the 
sensor may include a sensor body and a plurality of 
component sensors disposed on the sensor body, each 
component sensor being operative to measure the 
magnitude of a magnetic field component in a preselected 
local direction relative to the sensor body to provide a 
component sensor signal representing that particular 
component of the field prevailing at the sensor. The 
apparatus further includes calculation means for 
calculating the location of the sensor, the orientation 
of the sensor or both based upon the sensor signals and 
upon the unknown patterns of variation of the fields. 
The apparatus also includes feedback control means for 
adjusting the field generating means to alter the known 
pattern of variation of at least one of the fields 
responsive to the sensor signals, to the calculated 
location of the sensor or both so as to maintain the 
detected parameters of the altered field at the sensor 
within a preselected range. For example, where the 
field generating means includes a plurality of 
transmitters disposed adjacent the sensing volume, the 
feedback control means can be arranged to increase the 
strength of the field transmitted by a particular 
transmitter when the sensor is remote from such 
transmitter, and to decrease the strength of the field 
emitted by that transmitter when the sensor is close to 
the transmitter. Thus, the parameters of the field 
detected at the location of the sensor will always lie 
within a relatively narrow range of values. Because the 
transmitted field is altered in a known manner, the 
position and/or orientation of the sensor can still be 
calculated, simply by using the altered field as the 
basis for calculation. For example, where the feedback 
control means increases the current flowing to an 
electromagnetic coil, the parameters used in the 
calculation to represent the strength of the field 



WO 96/41119 

PCT/US96/08411 

emitted fro, that particular coil are adjusts 

accordingly . J ustea 

Apparatus according to this aspect of th. 
present invention can maintain the detected fi^d 
parameters at the sensor location, such as the magnitude 
of the field at such location, within a relatively 
narrow range even where the field strength varies with I 
power of distance greater than the first power as, ^or 
example, where the field strength varies with the th rd 
» Power of distance fro* the coil. The system ^J^* 
can provxde good signal to noise ratio, even with a ve" 
small sensor. Moreover, there is no need for any 
particular coil configuration to produce a linear or 
quasi-linear fie!d. For example, in a ZiZ 

disclosed at various locations beneath a patient 
receiving bed, or on one side of the bed. 

Although apparatus according to this aspect of 

- additional advantages when used with magnetoresist^ 
sensors or other sensors which can lose accuracy „h» 
exposed to very high Magnetic fi el ds. Because the fie" 
is maintained within a relatively narrow range, the 
sensors are not exposed to fields which would ^ 

Se f T7" ln 8 tUrthSr Varlant ° f invention' 
the feedback control system is actuated to maintain the 
parameters of the field at each sensor within ^ 
narrowest attainable range of values. Thus ^e 

30 se edba °\° 0ntr01 » aintai - the field parameters at the 
sensor at substantially constant v al ues. Here again! 

an*" "d £ ' eldS "~ altered in 3 ~ 

and the deposition of the sensor can still be 

35 

provide, fU ? her aSPSCt ° f thS Present invention 

adapted f - 1Ual " g " ^ ~ transmitter 

adapted for mounting to the object to be located and 

also Eluding a plurality of sensors disposed adjacent 
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a sensing volume. The apparatus according to this 
aspect of the invention includes calculation means for 
determining a calculated location of the transmitter 
based upon signals from the sensors and upon the known 
5 pattern of variation of the field at the transmitter. 
Here again, feedback control means are provided for 
adjusting the field generating means to alter the known 
pattern of variation of field sent by the at least one 
transmitter in response to the sensor signal, to the 

10 calculated location of the transmitter or both so as to 
maintain the parameters of the field from the at least 
one transmitter prevailing at each sensor within a 
preselected range. This arrangement is similar to the 
arrangement discussed above, except that the transmitter 

15 or transmitters are disposed on the object to be located 
whereas the sensors are located in the fixed frame of 
reference of the sensing volume. Here again, when the 
transmitter is close to the sensor, or where the 
orientations of the transmitter and sensor favor strong 

20 coupling between the transmitted field and the sensor, 
the strength of the field is reduced. If the 

transmitter is relatively far from the sensor, or in a 
unfavorable orientation, the control means increases the 
strength of the transmitted field. This arrangement 

25 provides advantages similar to those discussed above. 

A further aspect of the invention provides 
apparatus for sensing the disposition of an object — its 
position, orientation or both — in a frame of reference. 
Apparatus according to this aspect of the invention 

30 includes transmitter means including at least one 
transmitter for providing a field and sensor means 
including at least one sensor for sensing one or more 
parameters of a field prevailing at each such sensor and 
providing one or more sensor signals indicative of such 

35 parameters. The transmitter means and said sensor means 
cooperatively define a plurality of transmitter-sensor 
pairs, each including one transmitter and one sensor as 
elements of the pair, one element of each such pair 
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being disposed on the objeot and the other element or 
each such pair being disposea „ a ^ « 

said frame of reference. Typicallv , . 

e^ent o f each tra nOT itter-se„ Z ^ ~ 
5 different position or orientation *k ^ ^posed at a 
eie^ent of the o tter X ^ "'^ 
incxudes caicuiation ZZ £ 

rt: r a ob ^ b - -™ ™ 

signals feedback control means f«>- ^- . 
■0 fitting _ jTJTELJt 

P-vide T^Ltelinr ^™ 

» mentation wherein the ^ene^ 

means xs altered Vw 



20 



25 



30 



35 



aXtered by fe edbacK control based upon the 

L location 

™te rs o £ each t ieid to'-T^ected 

Within a nrocoi «~+.~,a bensor 



- * ^wiiuiux cased upon the 

£ zr ~- - — — — 



»u ik ,uj. atea location 

sensor of the t-ya«=™^4-*. 

pe . J 6 trans »itter, to maintain the 
rs of each field 

within a preselected range 
x>r M . * • 0thSr ° bjeCtS ' fe ^ures and advantages of the 

d ~ Wi11 bG ^ **- th 

forth bel T> ° f ^ -^diments set 

Lllgs ° W ' ^ C ° njUnCti0n ^ accompanying 

BRIEF DEfiOFTPTTOW .™ E DBMgT1MC , 

toinf- 1 ^ 9 diagraMnat i° Prospective view 

depicting portions of apparatus in accordance with Z 
embodiment of the invention. lth ° ne 

Figure 2 is a fragmentary, diam- amma -n 
prospective view depicting a portion of J™*™ 1 * 
illustrated in Fig. l ° f *** W^tus 

depictina V*T ' ** & fUnCtio »*l block diagram 

ZTTJT portions of the apparatus depictJ in 

steps in ^hU ^ 
-ention. Pigures 5 J d * ^^^0 
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perspective view depicting portions of apparatus in 
accordance with a further embodiment of the invention. 
BEST MODE OF CARRYING OUT INVENTION 

Apparatus in accordance with one embodiment of 
5 the present invention includes three generally helical 
transmitter coils 10 disposed in a common plane. Coils 
10 are mounted in fixed position in the frame of 
reference of patient receiving bed 12. That frame of 
reference is denoted by a Cartesian coordinate system X, 
10 Y f 2 as shown in Fig. 1. A patient P may be positioned 
on the patient receiving bed. The axes 11 of the coils 
are parallel to one another. A patient receiving bed 12 
extends just above the plane of coil 10. The apparatus 
further includes an object or probe 14. The probe is 
adapted to be inserted into a medical instrument such as 
a catheter 16 and positioned at a desired location in 
the catheter, such as at the distal tip of the catheter 
or at another location along the length of the catheter. 
Probe 14 has mounted thereon a sensor 18. Sensor 18 
includes three component sensors 20, 2 2 and 24 adapted 
to sense components of the magnetic field in mutually 
orthogonal local directions X f , Y ■ and Z'. That is, 
component sensor 20 is sensitive to magnetic fields 
directed in direction X', but largely insensitive to 
25 fields in directions Y 1 and Z«, whereas sensor 22 is 
sensitive only to fields in direction Y f and sensor 24 
is sensitive to fields in direction Z 1 . These sensors 
are adapted to provide separate sensor signals 
representing the separate components. Sensor 18 may be 
30 a solid state sensor of the type described in the 
aforementioned International Patent Publication 
W095/09562, the disclosure of which is hereby 
incorporated by reference herein. As further described 
therein, each of the component sensors may include a 
generally planar magnetically sensitive film, such as a 
magneto-resistive film or a Hall effect sensing film. 
Each such film may be sensitive to fields directed 
normal to the plane of the film. Alternatively, sensor 
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18 may include an array of miniature coils, the axes of 
the coxis being oriented orthogonal tQ ^ anothe ^ 

Although these represent the preferred sensors 
essentially any other magnetically sensitive device may 
be employed as, for example, magneto-optical sensors and 

riux gate magnetometers. 

component sensors 20, 22 and 24 are connected 
through a cable 26 , with separate leads for each sensor, 
to a command unit 28. Command unit 28 (Fig. 3 ) includes 
10 an input amplification and analog to digital ("a/D") 
conversion section 30 adapted to receive the individual 
signals from component sensors 20, 22 and 24 of sensor 
18, amplify the same and convert the same to digital 
form. The amplification and A/D conversion unit 30 may 
also include other conventional signal processing 
devices such as analog or digital band pass filtering 
and noise rejection devices and signal averagers. 
Command unit 28 further includes a computation unit 32. 

20 TT^ 32 ^ ^ implemen ted as a P-grammed 

general purpose computer. As further discussed below 
the position computation unit is arranged to compute the 
disposition of sensor i 8f and hence the disposition Qf 
the object or probe 14 at the catheter tip from the 
sensor signals. As used in this disclosure, the term 

25 "disposition" of an element refers to the position of 
the element, the orientation of the element or both 
Thus the computation unit is arranged to calculate the 
position of sensor 18, the orientation of the sensor 

30 un;t P 2 T rab1 "' b ° th POSiti ° n ^ command 
unxt 28 may be linKed to a dispiay device (not shown) 

for providing a human intelli gible representation of the 

t^elll^r ^ Pr ° be ^ ° bjeCt SUCh *«» 

intelligible representation may be provided either as 

35 Z^IT- inf ° rmati0n Presen ^ the position and/or 
orientation of object 14 in the X, Y, « coordinate 
system or, preferably, as a pictorial representation of 
the object and of the associated catheter superposed on 
a pictorial representation of the patient 
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Command unit 28 further includes a control 
unit 34. Control unit 34 is linked by output lines 36, 
38 and 4 0 to three separate coil drivers 42 , 44 and 46. 
Each coil driver is linked to a separate one of 
5 transmitter coils 10. Each coil driver is adapted to 
send a direct current through the associated transmit 
coil 10. Each coil driver is arranged to control the 
amplitude of such current, and to turn the current on or 
off, in response to control signals received from 
10 control unit 34. The control unit is arranged to signal 
the coil drivers to provide currents to their respective 
transmit coils in alternating sequence, so that coil 10a 
receives current while coils lob and 10c are inactive; 
coil 10b receives current while coils 10a and 10c are 
inactive and coil 10c receives current while coils 10a 
and 10b are inactive. The control unit receives data 
from amplification and conversion section 30 and, as 
further discussed below, actuates the coil drivers to 
vary the amplitude of the current to each coil. The 
control unit may include conventional interface devices 
such as digital to analog converters or bus interface 
units so that the output of the control unit is 
compatible with the control input of each coil driver. 
Also, although the control unit is illustrated 
separately from the other logical units of command unit 
28, it should be appreciated that the control unit may 
share physical elements of the command unit and other 
elements. For example, where the command unit 

incorporates a general purpose computer, the processor 
of the computer may serve both as an element of the 
position computation unit and as an element of the 
control unit, executing functions appropriate to the 
different units at different times. 

In a method according to one aspect of the 
invention, an catheter 16 is advanced into the body of a 
patient P. The probe 14 with sensor 18 thereon is 
disposed at the tip of the catheter. The catheter tip 
is disposed at an unknown location somewhere above the 
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Plane of coils io. control unit 34 actuates the coil 
drivers in sequence, using an initial or default value 
for the current amplitude to be provided to each coil 
10. Amplification and conversion unit 30 samples the 
signal from each of the component sensors 20, 22 and 24 

f lovT tr Se i eCted time ^ be * inn - g of current 

t mi zr :r coii - For exampie ' at * 

time after the beginning of current flow through 
transmit coil 10 «, unit 30 takes a sample of the signal 
from each of component sensors 20, 22 and 24, and 
converts the same to digital format. command unit 28 

!ndlld T^" 9 t0tal fiSld baSed these 

individual signals. The total field magnitude is: 

W=Vi*Ar +(* 22 s 22 ) 2 +(* M s 24 ) 2 

«a,H T" m '' Bl ° a ^ ^ magnitude of the magnetic 
field vector at sensor 18 at the time that coil i 0a is 

actuated; 

Strmo ^ , K2 ° ±S 3 SSnSltivlt y faot ° r relating the signal 
strength from sensor 20 to the magnetic field component 
along axis X' ; 

S 20 is the signal strength from sensor 20 
during the actuation; and 

K 22 , s 22 and K 24 and s 24 are similar sensitivity 
constants and signal strengths for the other sensors 22 
and 24. 

in like manner, the system actuates coils lob 

HereT ln °T T ' USiDg the ™t strength. 

Here agam, the system computes the magnitude of the 
total f.eia vector prevailing at the sensor during 
actuation of coil lob and independently computes th! 
magmtude of the total field vector prevailing at the 
sensor during actuation of coil ioc. 

As shown in Fig. 4 , after detecting the 
component signals and calculating the total field 
magnitude prevailing at the sensor during actuation of 
each coil, the control unit determines whether or not 
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all of the field magnitudes are within a preselected 
range of magnitudes. This preselected range is chosen 
to lie within the optimum operating range of sensor 18. 
Thus f the minimum field magnitude is selected to lie 
well above the noise threshold of the system and above 
the minimum sensitivity level of the sensor, where the 
maximum field level is selected to lie well below the 
maximum limit of linearity of the sensor, and well below 
the maximum field which the sensor can tolerate without 
loss of accuracy. For a typical magneto-resistive 
sensor, which is most accurate and repeatable when used 
with fields of less than about 4 Gauss, the preselected 
range of field magnitudes may be from about 1.0 to about 
2.5 Gauss. With typical Hall effect sensors, which are 
most accurate when used with fields over about 3 0 Gauss, 
the preselected range will be above about 3 0 Gauss. If 
the three field vector magnitudes observed during 
actuation of the three units 10a, 10b and 10c all lie 
within the preselected range, than the system proceeds 
to calculate the position and orientation of sensor 18, 
and hence the position and orientation of probe 14 and 
the associated tip 16 of the catheter using conventional 
position finding algorithms. For example, the 

mathematical methods disclosed in United states Patent 
4,710,708 for finding positions using a multiple 
transmitting or receiving stations and a multi-axis 
sensor can be employed. The disclosure of said '708 
patent is hereby incorporated by reference herein. 
Briefly, the magnitude of the fields in each of the 
local or sensor directions X', Y', z' represented by 
each of the component sensor signals from each of the 
component sensors 20, 22 and 2 4 is a function of the 
overall strength of the field from the coil (also 
referred to as the magnetic dipole moment of the coil) , 
the distance from the particular coil to the sensor and 
the sensor rotation angles, i.e., the angles between the 
local directions X», Y 1 and Z 1 and the coil frame of 
reference directions X, Y and Z. When three component 
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sensor readings are collected during actuation cf three 
separate coils and equated with the component strength 
expressed as a Action of location for fields from a 
particular coil, they form a system of nine equations 

inr t r th unknoMns (the x - y - 2 10Cati ° n ° f «» — « 

and the three rotation angles, . The derivation of these 
equations is set forth in Appendix A. That system of 
equations can be solved by iterative methods such as the 
method of Marquardt or the method of Broyden for least- 

Z?' S SOlUtl ° n ° f " system of 

nonlinear equations. The command unit then provides 
output indicating the position and orientation of Z 
sensor and hence indicating the position and orientation 
of probe 14 and of the distal tip of the catheter. 

If one or more of the field magnitudes is out 
of the preselected range, the system does not calculate 
the position and orientation. Ins tead, control unit 34 
changes the field strength of the coil or coils 
associated with the out-of-range field strength. Fo r 
» sample where probe 14 and sensor is are relative^ 
close to con 10 a, the total field magnitude detected 
when con 10a is actuated will be above the preselected 
range. control unit 34 therefore win oommana eoll 



» actuation cycle. Conversely, if the probe and sensor 

sense"d f £r0 " COil 10C ' " eld 

sensed during actuation of coil loc with the default of 

current value win be below the preselected range. 

30 to " there£0re »«l -struct coil driver 46 

to increase the current to coil ioc on the next 
actuation cycle. The coil drivers may be arranged to 
vary the overall field strength or dipole moment of each 
coil stepwise, as by varying the current stepwise. Each 
increase or decrease commanded by control unit 34 may be 
one step. Alternatively, the control unit can calculate 
an increase or decrease proportional to the degree to 
Which the total field magnitude deviates from a target 
value within the preselected range. Thus, a relative^ 
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large change can be made when the field magnitude is far 
outside of the range, whereas a smaller change can be 
employed when the field magnitude is close to the range 
or within the range. The correction process continues 
5 until all field strengths are within the preselected 
range, whereupon the system calculates position and 
orientation. After the system finds values of coil 
currents which produce field magnitudes within the 
range, subsequent actuation cycles use these currents. 

10 During operation, as the physician utilizes catheter 16, 
the position of the catheter tip, and hence the position 
of probe 14 and sensor 18 may change. Any such change 
may bring one or more of the field magnitudes outside of 
the preselected range, whereupon the system will 

15 readjust the currents to the coils once again. 

As the feedback control unit readjusts the 
currents to the coils, the altered values of current are 
translated into new values for field strengths from the 
individual coils which are used in the aforementioned 

20 position determining equations. In this manner, the 
system assures that the sensor is always exposed to a 
field having a magnitude within the preselected range 
wherever the sensor is placed within a sensing volume 50 
extending over a preselected region above the plane of 

25 coils 10. The exact size of sensing volume 50 will 
depend upon the breadth of the preselected field 
magnitude range and the dynamic range of coil drivers 
42, 44 and 46, i.e., the degree to which the coil 
drivers can vary the currents. The size of sensing 

30 volume 50 within which the sensor will receive fields 
within the preselected field magnitude range from all 
coils will also depend upon the positions of the coils. 
However, for a typical system having three coils spaced 
at vertices of an equilateral triangle with sides about 

35 40 cm long, the sensing volume includes a region 
extending upwardly about 60 cm from the plane of the 
coils. At the plane of the coils, the sensing volume 
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extends about 20 cm beyond the equilateral triangle 
bounded by the coils. ngxe 

In the description set forth above, only one 
ensor » med to determine the position ana 
5 locatxon of only one object. However, plural objects 
and plural sensors may be utilized. As described in "he 
aforementioned International Publication WO95/09562 a 



Plurality of fiducial markers may be utilized. Each 
fiducial marker may include a sensor body 52 a ta.53 
10 and a sensor 54 of the same tvno 

™™ 4- ^ ^ yPS as sensor 18 releasably 

connected to the _. y 



co„„ ected to ^ sensor bodY . The an ;;;; ic ™^ 

°™° n " nlt 30 - "™ected to each of th"e 
•ddxtion.1 sensors in the same manner aS to sensor is 
The t ag ot eaoh sensor may ^ afflxM -or is 

« Patxent. ana an image of the patient may be acquit 
usxng any waging method which will show ^ J. 

s„ 0 „ t TV f patient ' s ^ ana - hioh » l " 

show the tags 53 as, for example. magnetic resonance 
xmagxng (HEI) , X - ray , CAT scanning or the li*e. 
the tags are constructed to facilitate imaging in tta 
desxred im a gi„ g method. where x- ray imeging is used, 
the tag may be radioopaque; where MR1 Paging is used 

rtsonan 9 "* ^ 3 "^i* 1 " hl = h — -*-tic 

* bodllH reSP ° nSe readUy disti "^»°"* ««- that of 
bodily tissue, ordinarily, the image is acquired prior 
to use of the catheter and probe discussed above. Z 

steT S J t 4 IT 31 * ^ ^ PreMnt in «» 

step. After the waging step, the sensor bodies 54 and 

sensors of the fiducial makers are attached to Se 
30 patient's body at the positions of the tags. whence 
catheter is employed, sensor signals from sensor la on 
the probe in the catheter are acquired concurrently with 
signals from the sensor 52 on th. fl*, • , 

^ on rne fiducial markers. The 
system determines the location and orientation nr 
35 5a ar ,ri k ^ . orientation of sensor 

fields f MrkerS 52 usin * -9netic 

fields from cons 10 i„ the same manner as the system 
determines the position and orientation of sensor J and 
Probe 14 . The acquired position and orientation of the 
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fiducial markers 52 can be used to register the 
previously acquired image of the patient with the 
position and orientation data for sensor 18 and probe 
14. For example, where the position and orientation 
5 data for the sensor 18 and probe 14 are displayed as a 
pictorial image of the catheter tip 16 on a screen, the 
position and orientation of fiducial markers, as 
determined by the magnetic field locating method can 
also be displayed. The previously acquired image data 

10 may be displayed as a pictorial representation of the 
patient 1 s structures and that representation will also 
include a picture of the tags 53. The previously 
acquired images of the tags are brought into 

registration with the images of the fiducial markers 

15 derived from magnetic location by transforming one or 
the other image until the picture of the tags in the 
previously acquired image overlies the representation of 
the fiducial marker in the magnetically acquired image. 
Such registration can be achieved by manually adjusting 

20 inputs to an image display system while visually 
matching the displayed images, as described in the 
aforementioned "562 International Publication, or else 
by automatic calculation of the transformation and 
rotation parameters necessary to bring the images into 

25 registration. Such parameters can be calculated by 
equating the positions of three or more fiducial marker 
sensors in the magnetic location frame of reference (the 
X,Y,Z coordinates of Fig. 1) with the positions of the 
tags 53 in the image frame of reference as modified by 

30 a matrix of unknown rotations a further matrix of 
unknown transformations. As further described in 
Appendix B, the resulting matrix equation yields an 
overdetermined system of nonlinear equations which 
include the rotation angles and transformation 

35 distances. When the images are properly registered, the 
position of the catheter tip is displayed in true 
location and orientation relative to the patient's 
internal structures. In similar fashion, a greater 
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number of filial markers. or ^ 

instrument may be located. medical 
m a system using multiple sensors, such as 
sensors IB and 54, the coils may be operated in separate 

maantt T ^ Pr ° Viae 6a ° h SenS ° r " ith *i— ""tog 
magnitudes in the nreqpionf^ ng 
sensor w Preselected range appropriate for the 

sensor. For example, control unif 

. . ' u,J " tro -i- unit 34 may actuate the 

oil drivers 42 , 4, and «, M hence transmitting ^ 

.0 tt fi t , SeC ° nd ' thlrd "* fOUrth °uri„g 
and *" St . Cy ° le ' the «*- sensor i, are acquire! 

r g L S eT al E S frOT ~ S " ° f ^ 
are ignored. During the first +.u 

<c ^- ^ irst c ycle, the coil currentlv 

is adjusted to yield field magnitudes at sensors 
within the presented range as discussed above. D urin"g 

first fiducial marker are acquired, the signals from 
sensor is and the other fiducial marker sensors 54b, Z 
are ignored and the coil currents are adjusted to 
provide a field magnitude at sensor 54a within its 

rtr a ran9e - in the thira — - - » 

sensor 54b are acquired whereas the other signals are 
ignore; in the fourth cycle, the s i„r,= l. * 

, . .. J"", tne signals from sensor 54c 

are acquired and others are ignored. The cycles may be 
sequential, with the entirety of a first cyl e betoq 
performed before the entirety of a second cycle. The 
cycles may also be interspersed with one enoLer Zl 

T7~ oLl ariver 42 can he ~™ - p-i: 

a current appropriate to yield the proper field at the 
first sensor 18 and then can be actuated to provide the 
appropriate fields at the fiduol.i , 
54b r„ ■ fiducial marker sensors 54a, 

54b and 54c m succession, followed by a similar 
multiple actuation of the other coil. t-v sl " llar 

raier colls. Likewise, where 
more sensors are employed, there may be as many separa" 
cycles of operation as sensors. alternatively, if SO me 

I TT7 ^ knOKn ^ ^ Cl ° Se t0 ~r, " 

data from o^T """" ^ 

less thai r T"*' nUBter ° f »- 

less than the number of sensors. Thus, the currents 
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used on a given cycle may be adjusted to provide fields 
within the desired range of magnitude for one sensor. 
Signals from other sensor can be acquired during this 
cycle, although such other sensor may encounter a field 
5 outside the preselected range. In some cases, the coil 
currents are fixed for some cycles and adjusted for 
others. For example, where numerous fiducial markers 
are employed, it may be possible to achieve acceptable 
registration accuracy using fixed coil currents to 
10 provide the fields when data is acquired from the 
fiducial markers, whereas the coil currents may be 
adjusted in the manner discussed above when data is 
acquired from the sensor in the active device or 
catheter . 

15 Certain sensors tend to lose accuracy when 

exposed to magnetic fields above a predetermined 
maximum. For example, certain magneto-resistive sensors 
temporarily lose accuracy if exposed to magnetic fields 
above about 4 Gauss. Where such sensors are employed in 

20 a multi-sensor system, the system should have 
appropriate provisions to avoid exposing sensors to 
excessive fields. For example, if sensor 54a is 
disposed close to coil 10b, whereas sensor 18 is 
disposed remote from coil 10b, sensor 54a may be exposed 

25 to extremely high field levels if the system adjusts the 
current in coil 10b to the level required to generate 
the appropriate field at sensor 18c To avoid this 
condition, the data from all sensors can be acquired 
during all of the cycles, and control unit 34 can be 

30 arranged to increase coil currents in a progressive 
manner over several cycles when an increase is required. 
The data from sensors which are not used for position 
monitoring on a particular cycle can be used to inhibit 
further increases if the field at the nominally unused 

35 sensor is approaching dangerous levels. Thus, if the 
system is in the process of progressively increasing the 
coil current in coil 10b to provide an adequate field 
level at sensor 18, the system may terminate such 
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ZZTIV' T field nagn2i ° tude at ™ — « «• 

reaches toe * CY=le aSS ° Ciatea «ith first sensor is 

If th!l ™: X1Mn allowed at the second sensor. 

If this condition occurs while the field level at fi^ 
> sensor as is still belo „ ^ preselected 

to calculate position and orientation based on 
out-of-range sensor signal, or both, alternately if 
the accuracy of the sensor can be restore * k , 
» ft- the nominally unused sensors ca„ t -ed to 
initiate restoration. Po ^ , ed t0 

magnetoresistive sensors use s b a . eXample ' ce "ain 
If exnosad i-„ biasing magnetic field. 

rLt 7 eXcesslve "eWs. such sensors can be 

15 2 r fc reSt0red to — cy after the excess fieW is 
.5 removed by adjusting a bias magnetic field annUed 

~- ™" *™ -V he arrant t 

trigger the reset process ~ 

process for one sensor if it ( = 
exposed to excess field dur™ » "it is 

another sensor. * ' CyCle ass ° ci ated with 

20 In the foregoing discussion, the coils =-„ 

.riven alternately in . time-multiplexed arranglent 
F-guency domain multiplexing can also be empWd' 
Thus, each coil can be actuated with a driving™ 

» Tv di b ff : rent carrier fre ~- ™ — *™ - 
zi:; rim tci r ntinuous - The 

Afferent carrL S^LTST cT^ * ^ 

» ana log or digital f ilter'ing^hX. xh" 

the carrier frequency of coil i 0a can be used to 
r :? 6 t0tal " eld attributable to cotl 

ontrol unit 34 to adjust the current to coil i 0a The 

» ™ Lingte 001 " ° 0ilS ta ^ 

ca^nr^Tre h u:nc C i°r entS ^'fro ^ 
^iexing approach can le ex d^Ho ^ 



rr " 1 * €W1 De extended to multiple 
sensors. For example, each coil can be drivel at I 
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multiple carrier frequencies, equal to the number of 
sensors, the carrier frequencies for each coil being 
different from one another and different from all of the 
other carrier frequencies for all of the other coils. 
The components in the signals from the first sensor 
varying at the respective first carrier frequencies can 
be detected and the strength of the coil currents at the 
respective first carrier frequencies can be adjusted 
accordingly, so that the field magnitudes at the first 
carrier frequencies prevailing at sensor 18 will be 
within the desired range. In this stage of the process, 
the components varying at the other carrier frequencies 
are ignored. The process is reversed with respect the 
fiducial marker sensor 54a. This approach can be 
15 extended to greater numbers of sensors. Also, the 
frequency multiplexing and time multiplexing approaches 
can be combined. Thus, each coil can be driven at only 
one carrier frequency. During a first cycle, first 
sensor 18 is actuated and the components varying at all 
of the carrier frequencies are monitored. The coil 
currents are adjusted accordingly to bring the fields at 
each carrier frequency into the desired range. During a 
second cycle, the second sensor 54a (the first fiducial 
marker sensor) is used and the coil currents are 
adjusted in the same manner, and so on with the other 
fiducial marker sensors. 

In a further variant, the coil currents, and 
hence the strengths of the fields from the individual 
transmitters, can be adjusted to bring the field 
components detected by each individual component sensor 
to within a preselected range. In such a system, each 
coil is adjusted separately with respect to each 
component sensor. In a first cycle, the current to 
transmit coil 10a is adjusted to bring the individual 
sensor signal from component sensor 20, representing the 
field component magnitude in the local X' direction, 
into a preselected range. During this cycle, the 
signals from the other component sensors 22 and 24 are 
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disregarded. In the next = y =l e , coil l0> t . rea(Jjusted 
L° to 9 the / ield C ° mPOnent ln the PTl- direct 

. r e on component sensor 22 int ° «» 

range. The same coil 10a then brings the field 
component in local direction into . preselected 

::ii:„ a ;: the sl9nai froB the prim ° - 

section 24 is monitored rrv,,*„ „ 

' Thls se quence of operations is 

censor a9aln ea ° h ° £ ° th - transmitt r 

COllS, This annrna^K « - 



10 coils 



coil.. Thxs approach can be extended to plural sensors 
each having plural component sensors. Here again the 
system keeps tracK of the coil current magnitudes 
to produce the fields since b y each component sensor! 

» derail ^"T V ^ — «* 

(overall f le ld strength) terms of the b1«,,h-„ 
equations. simultaneous 

Where the local direction associated with a 
particular component sensor is orthogonal, or nearly 
orthogonal, to the direction of the field produced by I 
z° particular coil at the san«n>- i» t . 

brino sensor, it may be impossible to 

bring the component in that local direction into the 
preselected range without either exceeding the tlr^ 

oTtroVa T ^ " P ™"* * total™ 

so strong as to impair one of the other sensors. In 

this instance, however, at least one of the other 
component sensors will receive a component having a 
magnitude in ^ presaleotad ^ ^ ™* = 

this vanant, the sig^l from M of 

sensors may be monitored during a cycle associated Zl 

Ll .T" C ~ nt ~" ~' — =™ 

applied to the coil may be limited to avoid exposing any 

other unused component sensor to an excessive Lll 
component in its sensing direction. 

In a further variant or 
» r elected range of magnitudes for ^ fXltp'on^ 
» any particular direction is narrowed to include only 

optiml PreSelSCtea P-rerably within the 

optimum range of accuracy of the particular component 
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sensor. The feedback control system thus operates to 
adjust the coil currents until the field component 
magnitude is at such single value. The position and 
orientation are calculated in the same way as discussed 
5 above. This variant has the advantage that nonlinearity 
in the component sensor response cannot affect the 
accuracy of the system. Provided that it is known that 
a particular reading from the sensor corresponds to the 
preselected value of field component magnitude, 
10 deviation from a linear relationship between component 
magnitude and sensor reading at other values of the 
component magnitude will not impair the accuracy of the 
system . 

Certain field sensors exhibit so-called "off- 
15 axis sensitivity". That is, the transfer function or 
relationship between field component magnitude along the 
sensitive axis of a particular component sensor and the 
reading from that component sensor varies when a strong 
field component orthogonal to such axis is present. 
20 Off-axis sensitivity can be corrected by using the 
readings from two component sensors to evaluate the 
magnitude of the field perpendicular to the sensitive 
axis of the third component sensor, and using that 
magnitude to determine a correction factor to be applied 
to the reading from the third component sensor. 

In the system described above, when the coils 
are driven with direct current, amplification and 
conversion unit 3 0 samples the data from each component 
sensor at a preselected delay time after the inception 
30 of the current. in a variant of this system, the 
command unit adjusts the delay time from inception of 
the current to sampling depending upon the current 
applied in a particular actuation. After inception of 
the current, the current gradually rises to its steady 
35 state value, and hence the overall field strength or 
dipole moment likewise rises gradually. The changing 
field induces eddy currents in electrically conductive 
materials which may be present in or near the sensing 
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fields' H-T" ° UrrentS lndUOe * d «"°"al -gnetic 
fieids, which can induce erroneous sensor readings. The 
eddy currents decay when the rate of change in the till 

5 ^ " eld aPPr ° aCheS ltS ««• va ue 

The delay trme should he long enough for the current to 

n« to an appreciable portion of its steady state 

va ue, and long enough for eddy currents to dissipate 

sufficients that they do not cause appreciable r^rs 

■o i^rr deiay tine is — — — 

™e =y« em illustrated schematically in Pi g S 

That „, the probe or object i 14 to be tracked is 
» egu lp ped with a transmitters lle incorporating three 
mature coils 122 and 124 ^ on * 

~ aXeS - *" *° «■ -«rol system 

abL with T <nCt Sh<>Wn) SiBUar to 
above „ lt h reference to Fig. 3. The fixed reference 

frame system has three undirectional sensors ll0 a, Zl 

and hoc mounted in £ixe d position relative to the 

com» r elVi " 9 • S ™ «• —ted in a 

common plane. The common plane of the sensors extends 
generally vertically on one side of the patient 

from^ 9 ted - 11,6 SSnSin9 V ° 1Une 150 «*«- ™"y 

ThlS SyStem ^ USed in essentially the 

same manner as the system discussed above. Here again 
the current delivered to each coil, and hence the dlpole 
moment of the transmitted fieid associated with Hch 
coxl „ adjusted so that the magnitude of the field 

33 Z lZZtZ eaoh seneor is wlthin a — - — 

Nmer ° US oombinati °n= and variations of the 

Z "i* f i r uss r above can te uti "- a ^: 

departing from the present invention. By way of 
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example, the transmitting coils 10 of Fig. l can be 
disposed in a vertical plane arrangement as illustrated 
for the sensors in Fig. 5. Also, the number of 
transmitters and sensors can be varied. For example, as 
disclosed in the aforementioned International 
Publication WO94/04938, the system may include a sensor 
having three mutually orthogonal sensing axes and a set 
of three mutually orthogonal transmitting coils having a 
common center in the fixed frame of reference. 
Alternatively, the system may include a single uniaxial 
transmitting coil on the object to be located and three 
sets of sensors on the affixed frame of reference, each 
such sensor set incorporating three receiving coils or 
component sensors adapted to detect fields in a mutually 
15 orthogonal direction. The reciprocal arrangement — 
with three sets of mutually orthogonal transmitting 
coils and a single uniaxial sensor on the object to be 
tracked — is also usable. In general, the sensors and 
transmitters should define a plurality of transmitter- 
20 receiver pairs, each including one element on the object 
and one element on the fixed frame of reference. 

In a further variant, the system adjusts the 
output of the transmitters in response to the calculated 
disposition of the object being tracked, rather than 
directly in response to the component signals or total 
field magnitude signal. Thus, the system can initially 
operate with the default of current values; derive an 
initial reading of the position and orientation of the 
object and then use that initially determined position 
and orientation to calculate the desired setting for 
each coil to achieve the desired field levels at the 
sensor. Such desired setting is selected to yield the 
field within the desired magnitude range at the object, 
assuming that the object has the position and 
orientation found in the initial reading. On the next 
cycle, the so-calculated coil currents are utilized and 
the process is repeated. In a variant of this approach, 
the system can store a lookup table listing the 
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TOriat. coi! currants for various combinations of 
object pos ltl on and orientation. Using the ^ ' 
determined position ana orientation! the system 
retrieves appropriate coil-current values from ^ 
5 lookup table for use in the next cycle. 

In the embodiments discussed above, the sensor 
- associated with a catheter. The same system ca Z 
used with other medica! instruments as, for examoif 
endoscopes and surgical instruments. ihe svstem 

T b r pplied to deteMine the -^^> - 

other than medical instruments. For exampie, it can be 
used xn tracking an input device for a computer 

of the f I 3 thMe ° ther Varlati <™ «* combinations 

of ^f. , -vention, the foregoing description 

of the preferred embodiments should be taken by way of 
mustratxon rather than by way of limitation of the 
invention as defined by the claims. 

IHDDSTPT y APPT.TI JIPTT T TT 

20 

The invention can be applied in the medical 
and surgical fields as we n as in w mecllc al 

industry. ln the health «» 
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Appendix A: calculation of Position and Orientation 

5 Provided that we know the physical 

configuration of the field generator (transmitter) 
which is fixed in position during operation, the 
magnetic field detected by a sensor is a function of the 
position and orientation of the sensor. In our system 
10 field generator coils are stimulated sequentially. The 
field sensed by a sensor (3 component sensors per probe) 
can be expressed in terms of position in x, y, z and 
orientation a, p, y (roll, pitch and yaw respectively) 
i.e.: 

15 

B[sensor] [coil] = f [sensor] [coil] (x, y, 

Z/ <*/ P/ Y) 

Where [sensor] designates a particular sensor 
and [coil] designates a particular transmitter coil. 
20 If the real field the sensor measures when 

[coil] is on is B 1 [sensor] [coil] , then theoretically, 

B» [sensor] [coil] = B[sensor] [coil] 

25 

i.e. 

B» [sensor] [coil] - f [sensor] [coil] (x, y, 
z, a, p, y) - 0.0 

30 

Since we have 3 sensors and 3 coils, the total 
equations are 9 with 6 unknowns (x, y, z for probe space 
location, a, p, y for its orientation) . By applying non- 
linear least square method, we can find unique x, y, z, 
35 a, p, y for the probe. 

Above shows the general ideal of the 
algorithm. In detail: 
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Assume the orthogonal x,Y,2 reference 
system (magnetic location 
coordinate) is described by matrix 



coordinate 



Cartesian 



e,= 



Gm em 

Bm em em 

Cm em em 



the probe's orthogonal system is 



10 



s P n e P n e,\% 
e p22 e,2i 
e P 3i e P 32 e, 33 



and, since the three sensors on probe mav no* 
be orthogonal to each other their nnn IT 
can described as: non-orthogonal axes 



15 



C«n e„n e n n 

£«21 e«22 e„23 



a transfer matrix TMir-ii i- • 

t?ft . 1 ^ «^J.x J-LiJCD] which xs to be used 

» m later on calculation can be obtained ,ron>: 

^mm = e^ n • e rt/J 

V/J €{1,2,3} 

25 

another matrix orthoovrilMi . 
ue-ori , - v l x i I JJ which is to be 

used also can be defined as: 

- raw e m e w V/,y e {1,2,3} 
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since we use roll(oc), pitch(P), yaw(y) to 
define probe orientation, orthojDV [ i ] [ j ] can be also 
described by: 



10 



20 



sin(P) sin(y) 
sin(p) cos(y) 



15 sin(p) sin(y) 
sin(P) cos(y) 



ortho_ 


_OV[l][l] 


= cos (a) 


cos (y) 


- sin(a) 


ortho_ 


_OV[l][2] 


= cos (a) 


sin(y) 


- sin (a) 


ortho_ 


_OV[l][3] 


= -sin(a) 


cos(p) 




ortho_ 


_0V[2][1] 


= -cos(P) 


sin(y) 




ortho_ 


_0V[2] [2] 


- cos(p) 


cos (y) 




ortho_ 


.OV[2][3] 


= sin(p) 






ortho_ 


_OV[3][l] 


= sin (a) 


cos (y) 


+ cos (a) 


ortho_ 


OV[3][2] 


= sin (a) 


sin(y) 


- cos (a) 


ortho_ 


OV[3] [3] 


= cos (a) i 


cos(P) 





orthogonal vector matrix can therefore be 
calculated by matrix multiplication of previous defined 
matrix T and orhto OV: 



25 



ov= T*ortho OV 



The theoretical magnetic field for an 
orthogonal system generated at sensor position pointing 

to ^direction can be expressed as: 



30 f[™inifKx 9 y,z,a 9 J3,y) 

(detail abbreviated) , where f is a known function, and 
includes a dipole moment term having magnitude 
proportional to the current flow in the particular coil. 
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may „ ot n ° n -° rth °9° nal "y correction (sensors 

may not be perpendicular to each other,, the magnetic 
field sensor measures should be- 

5 



Assume the real field S ensor detected when 

10 coil is on is B '[sensor][coil] 

/ unen : 

B'[sensor][coH] - E[sensor][coil] « 0.0 

therefore, the 9 equations to be solved for x 
W Y, z, a, 0, Y , ar e: 



20 



*I««ni<af3] - 1 /[«-f3GPX»,, A «^ r ) x ov-p]^,] = o.o 



25 
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i=1 =0 . 0 



B'isemorlftcoiiq-jr/lcommix^z,^?,/) x ov ~ l \i][sensor2\ = 0.0 
j=i 

5 



B'[sensor3][coil\\ - ^f[coil\-\[i]{x,y,z, a,fi, r ) x ov '\i][serjsor3] = 0.0 



/=1 



£T^or3][co//2]~^/[co//2][/](x J ^r,a,Ar) x ov ~\i][semor3] = 0.0 



10 

3 

5 , [^or3][w/73]-^/[co/73][/](x J ^r,a,^r) * ov~ l [i][sensor3] = 0.0 



15 



A well-known non- linear least squares equation 
solver can be applied to solve the above equations and 
find probe position x, y, z and orientation a, p, y. 
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Appendix b: Registration of Fiducial Markers 

Assumption: 

nni „- • . " e kn °" **" coordinates of three fiducial 

ZLT, 1M9e data <MR eto -> coordinate ~ 

.0 My3 , {<MX1 ' Myl ' M21> ' <M * 2 < «»>. (MX3, 

fid„= a i • * ^ °°° rainates of the same three 

tiducal points in the x v 7 j- 

fixed fra.e of r^erel'ce of Th " W " °* 

" apparatus: ref — * of the n ac,neti= locating 



20 



Py3, Pz3)} 



Pyl, Pel), (Px2 , Py 2 , Pz2) , (Px3/ 



rx3 v, ^ points {(x1 ' y1 ' 21) - < x2 < «), 

(X3, y3 «)} can be formed by {(al , bl , o1) 

Th ' ?? thr ° Ugh r0tati °" * - d transition * 

ine rotation matrix R is: 



25 



iter iizx" 

toy tfyy toy 

toz i2y z Rzz 



where 



30 
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Rxx = cos(a) x cos(» - an(a) x sinQ?) x sin{» 
Rxy - cos(a) x sin(^) + sin(a) x sin(/?) x co^y) 
Rxz = - sin(a) x cos(/?) 
/Jyx: = - cos(/?) x sin(j) 

= cos(/?) x cos(y) 
i?yz = sin(/?) 

.Rzx = sin(a) x cos(» + cos(a) x sinQS) x sin(» 
Rzy = sin(a) x sinOO - cos(» x sin(/?) x cos^) 
.Rzz = cos(a) x cos(/?) 

The translation matrix T is: 



x y z 
x y z 
x y z 



The relationship between image position and 
magnetic location apparatus frame of reference position 



is: 



10 



Mxl 


Myl 


Mzl' 




~Px\ 


Py\ 


Pzl 


Mxl 


Myl 


Mzl 


R + T = 


Pxl 


Py2 


Pz2 


Mx3 


My3 


Mz3 




Px3 


Py3 


Pz3 



or 



Mxl 


Myl 


Mzl 




'Pxl 


Pyl 


Pzl 




"o 


0 


0' 


Mxl 


Myl 


Mzl 


R + T- 


Pxl 


Pyl 


Pzl 




0 


0 


0 


Mx3 


My3 


Mz3 




Px3 


Py3 


Pz3 




0 


0 


0 



15 

This forms 9 single equations: 
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Mx\xRxx + My\xRxy + Mz\xRxz + x-P x \ = 0 
MxlxR y x + My\xRyy + Mz\xRyz + y-Pyl = 0 

MxlxRzx + MylxRzy + MzlxRzz + z-Pzl = 0 

Mx2xRxx + My2xRxy + Mz2xRxz + x-Px2 = 0 

Mx2xRyx + My2xRyy + Mz2xRy Z +y-Py2 = 0 

Mx2xRzx + My2xRzy + Mz2xRzz + z-Pz2 = 0 
Mx3xRxx + My3xRxy + Mz3xRx 2+x -p x3 = 0 

M x2xRyx + My3xRyy + Mz3xRyz+y-.py 3 = 0 
Mx3xRzx + My3xRzy + Mz2xRzz + z-Pz3=:0 

(*) 

with 6 unknowns: a, B v fnr *. • 

' P/ ' tor rotation angles, 
x, y, z for space translation. 

The registration procedure is to find rotation 
-gles a, B, y and space translation x, y , z by solving 
the above 9 equations. 



10 



15 



20 



25 



in pseudo-code, the registration procedure is: 

registration () 
{ 

Initialize data buffer; 

Assign fiducial reference points to image 
data (usually 3 points) ; 9 

initialize magnetic locating system; 

*iau= ial point""" 6 P ° SitiOTS ° £ ^ *— 

^ square) 0 rzz^:, solvin9 nethoa <nm - iin - r 

and find rotation angles a, p, y and 
space translation x, y, 2; ' 

Form rotation matrix R from a, p, y and 
translation matrix T from x, y, z; 
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The rotation matrix R and translation matrix T 
can be applied to subsequent positions obtained, 
including the positions of the of the sensor associated 
with the catheter. 



5 
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CLAIMS: 

1. Apparatus for determining position 
comprising: 

(a) field generating means (i 0 ) for 
producing a plurality of magnetic fields, each said 
field having parameters varying with location within a 
sensing volume (50) according to a known pattern of 
variation, said patterns of variation being different 
for different ones of said fields; 

(b) at least one sensor (18) adapted to 
detect one or more of said parameters prevailing at such 
sensor when said sensor is at an unknown 

said sensing volume and to provide one or more sensor 
signals representing said one or more parameters; 

(c) calculation means (32 ) for 
determining a calculated location of said sensor based 

upon said sensor signals and «= a -i^ t 

y aiS and said known patterns of 

variation; and 

(d) feedback control means (34) for 
adjusting said field generating means to alter said 
known pattern of variation of at least one of said 

fields responsive to sairi 

said sensor signals, said 

cloal.t11oa.tion or both so as to maintain said 
palters ot each such altered Held prevailing at said 
sensor within a preselected range. 

2. Apparatus for determining position 

comprising: 
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(a) a plurality of sensors (110) adapted 
to detect one or more parameters of a magnetic field 
prevailing at each such sensor and to provide one or 
more sensor signals representing said one or more 

5 parameters, said sensors being disposed adjacent a 
sensing volume (150) at different positions, different 
orientations or both; 

(b) field generating means including at 
least one magnetic field transmitter (118) movable to an 

10 unknown location within said sensing volume and adapted 
to provide at least one magnetic field in said sensing 
volume, so that the magnetic field from each transmitter 
varies in according to a known pattern of variation with 
location relative to such transmitter whereby said 

15 sensor signals will represent one or more parameters of 
the field or fields from said at least one transmitter; 

(c) calculation means (32) for 
determining a calculated location of said at least one 
transmitter based upon said sensor signals and said 

20 known patterns of variation; and 

(d) feedback control means (34) for 
adjusting said field generating means to alter said 
known pattern of variation of at least one of said 
fields responsive to said sensor signals, said 

25 calculated location or both so as to maintain said 
parameters of each such altered field prevailing at each 
said sensor within a preselected range. 
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3. 



Apparatus as claimed in claim j. Q r olain 
2 wherein said field generating means is operative to 
provide DC Mgnetic tl . 1(ta wltMn ^ 

4. Apparatus as claimed in claim 3 wherein 
' said field generating means includes a plurality of 
coils (l0 . 120 , i22 , 124) and means for inter „ ittently 
actuating each said coil ,34. 42, 44, 46, with . direot 
current, said feedback control means includes means for 
altering the magnitude of each said direct current, and 

1» wherein said calculation means includes means for 
sampling the signals from each sensor when a delay time 
has elapsed after commencement of each said intermittent 
actuation, said feedback control means further including 
»eans for varying the delay time for each actuation in 

IS accordance with the magnitude of the current applied in 
such actuation so as to reduce the delay time when 
lesser currents are employed and increase the delay time 
when greater currents are employed. 

5- Apparatus as claimed in claim l or claim 

» 2 wherein said field generating means is operative to 
provide Ac magnetic fields within said sensing volume. 

6. Apparatus as claimed in claim l or claim 
2 wherein each said sensor has maximum accuracy within a 
preselected optimum operational range for each said 

» parameter and said preselected range for each said 
parameter substantially corresponds to said optimum 
operational range* 

7. Apparatus as claimed in claim l or claim 
2 Wherein said preselected range for each said parameter 
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consists of a single value, and wherein said feedback 
control means is operative to maintain each said 
parameter of the field at each sensor substantially at 
said single value. 
5 8. Apparatus as claimed in claim 1 wherein 

said field generating means includes a plurality of 
transmitters (10) disposed adjacent said sensing volume. 

9. Apparatus as claimed in claim 2 or claim 
8 wherein each said transmitter is operative to radiate 

10 a field having magnitude decreasing with r n where r is 
distance from the transmitter and n is a number greater 
than 1. 

10. Apparatus as claimed in claim 9 wherein 
each said transmitter is operative to radiate a field 

15 having magnitude substantially equal to K/r 3 where K is 
a real number, and wherein said feedback control means 
is operative to adjust K. 

11. Apparatus as claimed in claim 1 further 
comprising a sensor body (18), said at least one sensor 

20 including a plurality of component sensors (2 0, 22, 24) 
disposed on said sensor body, each said component sensor 
being operative to measure the magnitude of a magnetic 
field component in a preselected local direction 
relative to the sensor body. 

25 12 • Apparatus as claimed in claim 11 wherein 

said field generating means includes a plurality of 
transmitters (10) disposed at spaced-apart locations 
adjacent said sensing volume, said feedback control 
means being operative to adjust each said transmitter to 
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transmit a plurality of field, so that for each said 
field, the magnitude of a component of such field in 
the local direction of one said component sensor will be 
within a preselected range. 

13. Apparatus as claimed in claim n wherein 
each said component sensor is a H all effect sensor, a 
xnagnetoresistive sensor, a magnetooptical sensor or a 
flux gate magnetometer . 

14. Apparatus as claimed in claim i or claim 
10 2 wherein each said sensor includes a coil and wherein 

said field generating means is arranged to provide Ac 
fields in said sensing volume. 

15. Apparatus as claimed in claim i wherein 
said field generating means includes a plurality of 

15 transmitters disposed at spaced-apart locations on one 
side of a reference plane and wherein said sensing 
volume extends to the other side of said reference 
plane. 

Apparatus as claimed in claim i S wherein 
said transfers are disposed substantially in a common 

plane. 

17. Apparatus as claimed in c l aim 2 wherein 
said plurality of sensors are disposed at spaced-apart 
locations on one side of a reference plane and wherein 
said sensing volume extends to the other side of said 
reference plane. 

18. Apparatus as claimed in claim 17 wherein 
said sensors are disposed substantially in a coiranon 
plane. 



25 



20 

plane. 
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19. Apparatus as claimed in claim 1 or claim 
2 wherein said at least one sensor includes a plurality 
of sensors, and wherein said sensors can be impaired by 
exposure to fields having components above a 
5 predetermined threshold magnitude, and wherein said 
feedback control means includes means for controlling 
said field generating means to assure that no sensor is 
exposed to a field component having magnitude above its 
threshold value. 

10 20 • Apparatus as claimed in claim 1 or claim 

2 wherein said at least one sensor includes a plurality 
of sensors and wherein said sensors can be impaired by 
exposure to fields having components above a 
predetermined threshold magnitude but can be restored 

15 after such exposure upon application of a reset routine, 
said feedback control means including means for 
initiating such reset routine in the event of such 
exposure . 

21. Apparatus for sensing the disposition of 
10 an object in a frame of reference comprising: 

(a) transmitter means (10) including at 
least one transmitter for providing a field; 

(b) sensor means (18) including at least 
one sensor for sensing one or more parameters of a field 

15 prevailing at each such sensor and providing one or more 
sensor signals indicative of such parameters, said 
transmitter means and said sensor means cooperatively 
defining a plurality of transmitter-sensor pairs, each 
including one transmitter and one sensor as elements of 



WO 96/41119 

PCT/US96/08411 

-42- 

the pair, one element of each such p air being disposed 
on the object and the other element of each such pair 
being disposed at a Known disposition in said frame of 
reference ; 

5 (c) calculation means {32 ) for 

determining the disposition of said object in said frame 
of reference based upon said sensor signals; and 

(d) feedback control means (34 ) for 
adjusting said transmitting raeans responsive to said 
- sensor signals to thereby maintain at least one of said 
sensor signals within a preselected range. 

22. a method of determining the disposition 
of an object in a frame of reference comprising the 
steps of: 

(a) P r ° v i<*ing one or more magnetic 
fields each having a known pattern of variation with 
respect to said frame of reference; 

(b) sensing one or more parameters of 
said one or more magnetic fields prevailing at said 

» object to provide one or more sensor signals indicative 
of said one or more parameters; 

(c) calculating the disposition of said 
object in said frame of reference based upon said sensor 
Signals and said known pattern or patterns of variation- 

25 and 

(«> altering said known pattern of 

variation of at least one of said ««o 

or said one or more magnetic 

fields responsive to ca^ 

to said sensor signals, said 

calculated disposition or both so a «, „ • ^ . 

tn so as to maintain said 
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parameters of each such altered field prevailing at said 
object within a preselected range, 

23. A method of determining the disposition 
of an object in a frame of reference comprising the 
5 steps of : 

(a) providing one or more magnetic 
fields each having a known pattern of variation in space 
with respect to object; 

(b) sensing one or more parameters of 
10 said one or more magnetic fields prevailing at one or 

more sensors having known disposition in said frame of 
reference to provide one or more sensor signals 
indicative of said one or more parameters; 

(c) calculating the disposition of said 
L5 object in said frame of reference based upon said sensor 

signals and said known pattern or patterns of variation 
of said one or more magnetic fields; and 

(d) altering said known pattern of 
variation of at least one of said one or more magnetic 

!0 fields responsive to said sensor signals, said 
calculated disposition or both so as to maintain said 
parameters of each such altered field prevailing at at 
least one of said sensors within a preselected range. 
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